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Crystalline structures, zone-center phonon modes, and the related dielectric response of the three 
low-pressure phases of Hf02 have been investigated in density-functional theory using ultrasoft pseu- 
dopotentials and a plane-wave basis. The structures of low-pressure HfC>2 polymorphs are carefully 
studied with both the local-density approximation (LDA) and the generalized gradient approxima- 
tion (GGA). The fully relaxed structures obtained with either exchange-correlation scheme agree 
reasonably well with experiment, although LDA yields better overall agreement. After calculating 
the Born effective charge tensors and the force-constant matrices by finite-difference methods, the 
lattice dielectric susceptibility tensors for the three HfC>2 phases are computed by decomposing the 
tensors into the contributions from individual infrared-active phonon modes. 

PACS numbers: 77.22.-d, 61.66.-f, 63.20.-e, 77.84.Bw 



Hafnia (Hf02) is technologically important because of 
its extraordinary high bulk modulus, high melting point, 
and high chemical stability, as well as its high neutron 
absorption cross section. H1O2 resembles its twin oxide, 
zirconia (Zr02), in many physical and chemical proper- 
ties. The resemblance is attributable to the structural 
similarity between the two oxides, which can in turn be 
explained by the chemical similarity of Hf and Zr, which 
have similar atomic and ionic radii (i.e., ionic radii for 
Hf 4 + and Zr 4 + of 0.78 and 0.79 A, respectively §) as 
a result of the so-called lanthanide contraction. Under 
ambient pressure, both oxides are monoclinic (to, space 
group P2i/c) at low temperature, and transform to a 
tetragonal structure (f, space group PA2/ rune) and then 
to a cubic structure (c, space group Fm3m) as the tem- 
perature increases, as illustrated in Fig. [l| 

High-K metal-oxide dielectrics have recently been the 
focus of substantial ongoing efforts directed toward find- 
ing a replacement for SiC>2 as the gate dielectric in com- 
plementary metal-oxide-semiconductor (CMOS) devices. 
Hf02, Z1O2 and their SiC>2 mixtures show promise for 
this purpose [0J^]. Thus, a systematic theoretical investi- 
gation of the structural and dielectric properties of these 
dielectrics, in both bulk and thin-film form, is clearly de- 
sirable. As a first step in this direction, we have, in a 
previous paper 0], investigated the bulk structures and 
lattice dielectric response of ZrC>2 polymorphs. We found 
that the dielectric responses vary dramatically with the 
crystal phase. Specifically, we found that the monoclinic 
phase has a strongly anisotropic lattice dielectric ten- 
sor and a rather small orientationally-averaged dielectric 
constant owing to the fact that the mode effective charges 
associated with the softest modes are relatively weak. 

This report presents the corresponding work on HTO2, 
providing the first thorough theoretical study of the 
structural, vibrational and lattice dielectric properties of 
the HfC>2 phases. Such properties are naturally expected 
to be similar to those of ZrC>2 in view of the chemical sim- 



ilarities mentioned above. We find that this is generally 
true, although we also find some significant quantitative 
differences in some of the calculated properties. 

The calculation of the lattice contributions to the 
static dielectric tensor eo entails the computations of 
the Born effective charge tensors Z* and the force- 
constant matrices 3>. The Z* tensors, defined via 
AP = (e/V) J2t z * ■ Au * i are obtained by finite dif- 
ferences of polarizations (P) as various sublattice dis- 
placements (iij) are imposed, with the electronic part 
of the polarizations computed using the Berry-phase ap- 
proach |5|,|6|]. Here V is the volume of the unit cell, 
e is the electron charge, and i labels the atom in the 
unit cell. We then calculate the force-constant ma- 



trix, $ 



a/3 



- dF? /du? ~ - AF™ /An? by calculating all 
the Hellmann-Feynman forces F° caused by making dis- 
placements Uj of each atom in each Cartesian direction in 
turn (Greek indices label the Cartesian coordinates). The 
resulting $ matrix is symmetrized to clean up numerical 
errors, the dynamical matrix = (M^Af.,) -1 / 2 is 
constructed, and the latter is then diagonalized to obtain 
the eigenvalues uo\ and eigenvectors 

The static dielectric tensor can be decomposed into a 
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FIG. 1. Structures of the three HTO2 phases. Small dark 
circles and larger open circles denote Hf and O atoms respec- 
tively. Hf-O bonds are only shown in m-Hf02- In t-Hf02, 
the arrows indicate the shift of oxygen pairs. 
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TABLE I. Calculated structural parameters for three HfC>2 
phases using both LDA and GGA. Lattice parameters a, b, c 
are in A, (3 is in degrees, and V (volume per formula) is in 
A . Internal coordinates x, y and z are dimensionless. 



TABLE II. LDA dynamical effective charges Z* for HfC>2 
phases. (Values in parentheses are GGA results.) 
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contribution arising from purely electronic screening, 
and the contributions of the IR-active phonon modes, 
according to M 
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Si/3 ^i,a/: 



1/2 

(Mo/Mi) ' & t \p are mode 
effective charges, e is the electron charge, Mq is a refer- 
ence mass that we take for convenience to be 1 amu, u>\ is 
the frequency of the A-th IR-active phonon mode, and V 
is the volume of the 3-atom, 6-atom, or 12-atom unit cell 
for cubic, tetragonal, or monoclinic cases, respectively. 
£i,A/3, the eigendisplacement of atom i in phonon mode 
A, is normalized according to J2ia €i,\'a = S xx' ■ 

The calculations are carried out within an ultrasoft 
pseudopotential [g implementation of density-functional 
theory with a plane- wave basis and a conjugate-gradient 
minimization algorithm. The crystal structures of HfC>2 
polymorphs are investigated in the local-density approx- 
imation (LDA) as parameterized by Ceperley and Alder 
H as well as in the generalized gradient approximation 
(GGA) using PBE parametrization |I(J. We find that 
LDA yields slightly better agreement with the experi- 





Hf 




Oi 


o 2 




Cubic 


5.85 




-2.93 


-2.93 


Tetrag 


onal 










x'x' 


5.84 




-3.53 


-2.31 


1 1 

y y 


5.84 




-2.31 


-3.53 


zz 


5.00 




-2.50 


-2.50 


Monoclinic 










XX 


5.56 ( 


5.57) 


-3.09 (-3.10) 


-2.48 


(-2.47) 


xy 


-0.47 ( 


-0.56) 


0.97 ( 0.90) 


0.20 


( 0.15) 


xz 


0.96 ( 


0.91) 


-0.58 (-0.53) 


-0.39 


(-0.36) 


yx 


-0.13 ( 


-0.02) 


1.37 ( 1.29) 


0.21 


(0.11) 


yy 


5.55 ( 


5.57) 


-2.73 (-2.70) 


-2.82 


(-2.87) 


yz 


0.14 ( 


0.07) 


-0.71 (-0.61) 


0.35 


( 0.40) 


zx 


0.21 ( 


0.27) 


-0.18 (-0.20) 


-0.07 


(-0.09) 


zy 


0.41 ( 


0.45) 


-0.61 (-0.51) 


0.43 


( 0.46) 


zz 


4.74 ( 


4.64) 


-2.24 (-2.16) 


-2.58 


(-2.52) 



mental structures, and we therefore suggest that our LDA 
results for the dielectric properties are more reliable. The 
4s and Ap semicore shells are included in the valence in 
the Hf pseudopotential, and an energy cutoff of 25 Ry 
is chosen. A 4x4x4 Monkhorst-Pack k-point mesh is 
found to provide sufficient precision in the calculations 
of total energies and forces, and a 4x4x20 k-point sam- 
ling is used for calculating the Berry-phase polarization 
Each atomic sublattice is displaced in turn along 
each Cartesian direction by ±0.2% in lattice units, the 
electronic polarization and Hellmann-Feynman forces are 
computed, and Z* and are then constructed by finite 
differences from the results. 

Tabulated in Table | are the relaxed structural parame- 
ters for the three Hf02 polymorphs, with the correspond- 
ing data for Zr02 listed in the last column for comparison 
. While several structural determinations for m-Hf02 
can be found in the literature jl 11 1, corresponding re- 
sults for the tetragonal and cubic phases are relatively 
sparse JL2[ . Nor has there been much theoretical work 
on hafnia; most important is the recent work of Ref. |l5| ] 
which agrees quite well with our results. For m-Hf02, 
the parameters given in Ref. JHJ were used as the start- 
ing point of our relaxation procedures, while for t- and 
c-Hf02 we started the relaxation from the zirconia exper- 
imental structures. It can readily be seen that both the 
LDA and GGA agree reasonably well with the previous 
work, but that the LDA yields a better overall agreement. 
Our total-energy calculations reproduce the correct ener- 
getic ordering of the phases (monoclinic then tetragonal 
then cubic) using either LDA or GGA. 

Our results for the dynamical effective charges are pre- 
sented in Table [nj The symmetry of c-Hf02 requires 
that Z* be isotropic on each atom. In t-Hf02, the shift- 
ing of oxygen atoms creates two different configurations 
for oxygen atoms (denoted Oi and O2) and introduces 
off-diagonal elements. Thus, it is more natural to refer 
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TABLE III. Theoretical (LDA and GGA) and experimen- 
tal (Ref. [jl5|) frequencies (in cm" 1 ) of Raman- active phonon 
modes in monoclinic Hf02- 
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to a reference frame x'-y'-z that is rotated 45° about 
the z axis from the original Cartesian frame. Z*(Oi i 2) 
become diagonal in this frame. In rn-Hf02, there are 
two non-equivalent oxygen sites (i.e., the 3-fold and 4- 
fold oxygens, labeled as Oi and O2 respectively). The 
crystal structure can then be regarded as composed of 
three kinds of atoms, namely, Zr, Oi, and O2, which all 
have equally low symmetry, and their resulting Z* ten- 
sors are neither diagonal nor symmetric. The presence of 
two non-equivalent oxygen atoms with very different en- 
vironments is reflected in the difference between the Born 
effective charge tensors for Oi and O2. The anomalously 
large Z* values indicate that there is a strong dynamic 
charge transfer along the Hf— O bond as the bond length 
varies, indicating a mixed ionic-covalent nature of the 
Hf— O bond. The resultant relatively delocalized distri- 
bution of the electronic charge is very similar to ZrC>2, 
and is quite common in partially covalent oxides. 

Since Hf02 is isomorphic to Zr02, the analysis of the 
phonon modes at T is the same for HfOa as for Zr0 2 @. 
Of 36 phonon modes predicted for m-Hf02, 18 modes 
(9A g + 9B g ) are Raman-active and 15 modes (8A U + 7B U ) 
are IR-active, the remaining three modes being the zero- 
frequency translational modes. There are three IR-active 
modes {A2 U and two E u ) and three Raman-active modes 
(Ai g , Big and E g ) for £-HfC>2. Only one IR-active mode 
(one T\ u triplet) is predicted for C-HIO2. 

The Raman spectra of m-Hf02 have been extensively 
measured experimentally p^-^7|] , but the situation is not 
entirely satisfactory |]l6f . Issues concerning the number of 
modes and the mode assignments still remain unresolved, 
partially because of sample impurities and the broadness 



TABLE IV. Frequencies uj\ (cm 1 ) and scalar mode effec- 
tive charges Z\ of IR-active phonon modes for Hf02 phases, 
where = Y. a Ztl- 
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and weakness of some observed features. Thus, our ab- 
initio theoretical calculation can play an important role 
in establishing the Raman assignments. Table III shows 
the frequencies of the A g and B g Raman-active modes 
as calculated in LDA and GGA, together with the ob- 
served frequencies from a polarized Raman measurement 
on a high-quality single crystal jll|. The agreement is 
generally excellent; the observed Raman shifts mostly 
fall comfortably in the LDA-GGA range. A later single- 
crystal (but unpolarized) Raman spectrum |l7) shows al- 
most identical mode frequencies. However, a few details 
about the Table deserve comment, (i) We omit the weak 



mode reported as A g at 268 cm in Ref. p5| because 



and it does not fit with 
We assign the 323 cm _1 
as A g . (iii) The feature 
fL5| is presumed to be a 



it is not confirmed in Ref. |17fl 
our theoretical assignments, (ii) 
mode observed in Refs. Jl5| , p7t 
observed at 872 cm -1 in Ref. [ 
two-phonon process and is omitted, (iv) A weak mode 
is observed at 773 cm -1 in Ref. [fl7|| ; since this is consis- 
tent with our highest-frequency B g mode, we assign it as 
such. 

The frequencies of the IR-active phonon modes for the 
three Hf02 phases are tabulated in Table [V, together 
with the scalar mode effective charges. It can be seen that 
the frequencies calculated in GGA are shifted to lower 
frequency by ~10-16% relative to the LDA ones, while 
the mode assignments coincide exactly. As indicated in 
Eq. (Q), the contribution of a given IR-active mode to the 
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static dielectric constant scales as 2* 2 /w^ (§, so that one 

or more low- frequency modes with large Z* 's are needed 
to yield a large dielectric constant. As can be seen from 
Table [V, however, the few softest modes (< 300 cm" 1 ) 
have relatively small Z*'s, while the more active infrared 
modes come in the intermediate range of the IR spectrum 
(350 -450 cm -1 ). The general pattern is very similar to 
the case of ZrC>2. 

The lattice contributions to the dielectric tensors are 
obtained by summing the second term of Eq. (Q) over all 
the IR-active modes. Using the LDA we find 
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(The corresponding matrix elements of e l ^ no in the GGA 
tend to be larger than the LDA results by ~18%.) When 
compared with Zr02, the off-diagonal elements of e™;*^ 
are roughly doubled, while the diagonal elements become 
smaller. Most surprisingly, the x-y components of 4ct*a 
become more than twice as large as for Zr02, while the 
z component decreases by ~28%. We find the isotropic 
e cubic to be 23.9, somewhat smaller than the value of 31.8 
for Zr0 2 §. 

A direct comparison of these dielectric tensors with 
experiment is not feasible since there are few experimen- 
tal measurements, especially on the cubic and tetrago- 
nal phases. On the other hand, most measurements of 
which we are aware have been carried out on thin films 
(presumed to be monoclinic), and the reported dielectric 
constants span a wide range of 16-45 P, |l~8| , |i~9[ . Assum- 
ing an isotropic ~ 5 M], we obtained orientationally 
averaged static dielectric constants of 29, 70, and 16 (18 
in GGA) for the cubic, tetragonal and monoclinic Hf02 
phases, respectively. Our results then agree reasonably 
well with the more recent results in Ref. JT^j (thin film 
- 1700 A) and Rcf. g (ultrathin film < 100 A) which 
report eo to be 16 and 20 respectively. The surprising 
high eo measured in other experiments could possibly be 
explained by the presence of i-Hf02, which is known to 
be a metastable phase and which might be stablized by 
film stress or grain-size effects @[l9||20) 



In summary, we have investigated here the Born effec- 
tive charge tensors, zone-centered phonons, and the lat- 
tice contributions to the static dielectric tensors for the 
three Hf02 phases. It is found that the cubic and tetrag- 
onal phases have much larger dielectric response than the 



monoclinic phase, with an even stronger anisotropy in t- 
Hf02 than in t-ZiO-2- The overall dielectric constants for 
c-Hf02 and m-Hf02 are found to become smaller, while 
i-Hf02 has a much greater dielectric constant, than the 
corresponding values in Zr02- Moreover, our Raman re- 
sults can be used in resolving the puzzles associated with 
the Raman spectrum for 77i-Hf02. 

This work has been supported by NSF Grant DMR- 
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